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Abstract 

Differential dijet cross sections in diffractive deep-inelastic scattering are measured with 
the HI detector at HERA using an integrated luminosity of 51.5 pb~^. The selected events 
are of the type ep — > eXY, where the system X contains at least two jets and is well 
separated in rapidity from the low mass proton dissociation system Y. The dijet data are 
compared with QCD predictions at next-to-leading order based on diffractive parton dis- 
tribution functions previously extracted from measurements of inclusive diffractive deep- 
inelastic scattering. The prediction describes the dijet data well at low and intermediate zjp 
(the fraction of the momentum of the diffractive exchange canied by the parton entering the 
hard interaction) where the gluon density is well determined from the inclusive diffractive 
data, supporting QCD factorisation. A new set of diffractive parton distribution functions 
is obtained through a simultaneous fit to the diffractive inclusive and dijet cross sections. 
This allows for a precise determination of both the diffractive quark and gluon distributions 
in the range 0.05 < zjp < 0.9. In particular, the precision on the gluon density at high 
momentum fractions is improved compared to previous extractions. 
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1 Introduction 



Hadron-hadron collisions proceed predominantly via soft interactions to which perturbative 
Quantum Chromodynamics (QCD) cannot be applied. In a sizeable fraction of these soft pro- 
cesses the colliding hadrons remain intact or merely dissociate to larger mass states with the 
same quantum numbers. These "diffractive processes" dominate the behaviour of the total cross 
section at high energy and are phenomenologically described by the exchange of the pomeron 
trajectory, which carries the quantum numbers of the vacuum. The parton composition of this 
diffractive exchange is, however, not well known. 

Processes of the type ep eXp have been studied in detail at HERA. These processes can 
be pictured as 7*p scattering, where the virtual photon interacts with a diffractive exchange and 
dissociates to produce a system X. In QCD a hard scattering coUinear factorisation theorem [1] 
predicts that the cross section for diffractive deep-inelastic ep scattering (DIS) factorises into 
a set of universal diffractive parton distribution functions (DPDFs) of the proton and process- 
dependent hard scattering coefficients. DPDFs have been determined through QCD fits to the 
measured cross sections of inclusive diffractive scattering at HERA [2-15]. 

If QCD factorisation is fulfilled, next-to-leading order (NLO) QCD calculations based on 
DPDFs such as those extracted in [5] should be able to predict the production rates of more 
exclusive diffractive processes such as dijet and open charm production. Previous measure- 
ments of such exclusive cross sections in DIS [16-23] support QCD factorisation since they 
can be reasonably well described using the DPDFs determined from inclusive diffractive scat- 
tering. Diffractive dijet and charm production proceed mainly via boson gluon fusion (BGF, 
depicted in figure H]) and are therefore mainly sensitive to the diffractive gluon density. It was 
recently shown [5] that inclusive diffractive scattering data do not constrain the diffractive gluon 
density well at high momentum fractions. Thus stringent tests of factorisation can only be per- 
formed at low momentum fractions. However, the gluon density at high momentum fractions is 
particularly relevant for the estimation of cross sections for several important processes at the 
LHC [24]. Measurements of diffractive dijet production can directly constrain the diffractive 
gluon density at high momentum fractions, extending the kinematic range of reliably deter- 
mined diffractive parton densities. 

In this paper, a new measurement of diffractive dijet cross sections in deep-inelastic scat- 
tering is presented, based on data collected with the HI detector at HERA in the years 1999 
and 2000. These are the first HERA diffractive DIS dijet data with Ep = 920 GeV. Jets are 
defined using the inclusive k-r algorithm [25]. The resulting dijet cross sections are compared 
to NLO QCD predictions based on DPDFs previously extracted [5] from inclusive diffractive 
ep scattering at HI . For the first time, a combined NLO QCD fit is performed to the differential 
dijet cross sections and the inclusive diffractive cross section data in order to determine a new 
set of DPDFs. 



2 Kinematics 

The dominant process leading to the production of dijets in diffractive DIS is depicted in 
figure [TJ The incoming proton of four-momentum P interacts with the positron of four-momen- 
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Figure 1: Leading order diagram for diffractive dijet production in DIS. 



turn k via the exchange of a virtual photon with four-momentum q. The DIS kinematic variables 
are defined as 

n2 _ 2 _ P ■ Q 

^ ^ ' 2P-g' ^ P-k' 

where is the photon virtuality, x is the longitudinal momentum fraction of the proton carried 
by the struck quark and y is the inelasticity of the process. These quantities are connected by 
the relation 

= xys, 

where s denotes the fixed ep centre-of-mass energy squared. 

The hadronic final state of the events is divided into two systems X and Y, separated by 
the largest gap in the rapidity distribution of the hadrons relative to the collision axis in the 'j^p 
centre of mass system. The diffractive scattering is described in terms of the variables 

t={P-pYf, X]p= ^ J^^^ , P = x/X]p, 

q ■ F 

with py representing the four-momentum of the system Y . Here t is the squared four-momentum 
transfer at the proton vertex and xjp is the fraction of the proton's longitudinal momentum 
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transferred to the system X. The fractional longitudinal momentum of the diffractive exchange 
carried by the parton which enters the hard interaction with four-momentum v is given by 

q ■ V 
q-{P- Py) 

3 Experimental Procedure 

3.1 HI Detector 

A detailed description of the HI detector can be found in [26-28]. Here, a brief account of 
the components most relevant to the present analysis is given. The origin of the HI coordinate 
system is the nominal ep interaction point. The direction of the proton beam defines the positive 
2-axis (forward direction). Transverse momenta are measured in the x-y plane. Polar (9) 
and azimuthal (0) angles are measured with respect to this reference system. The pseudorapidity 
is defined as 77 = — lntan(^/2). 

The ep interaction region is surrounded by a two-layered silicon strip detector [29] and two 
large concentric drift chambers, operated inside a 1.16 T solenoidal magnetic field. Charged 
particle momenta are measured in the pseudorapidity range —1.5 < r] < 1.5 with a reso- 
lution of a{pT)/pT = 0.005 pr/GeV © 0.015. The central tracking detectors also provide 
triggering information based on track segments measured in the r-0 plane of the central jet 
chambers and on the z position of the event vertex obtained from the double layers of two 
multi-wire proportional chambers. A finely segmented electromagnetic and hadronic liquid 
argon (LAr) calorimeter [30] covers the range —1.5 < rj < 3.4. The energy resolution is 
a/E = 0.11/ ^jE/GeY for electromagnetic showers and a/E = 0.50/ ^/E/ GeV for hadrons, 
as measured in test beams [31]. A lead/scintillating fibre calorimeter (SPACAL) [28] covers the 
backward region —4 < 77 < — 1.4. Its main purpose is the detection of scattered positrons. 

The luminosity is measured via the Bethe-Heitler Bremsstrahlung process ep ep'j, the 
final state photon being detected in a crystal calorimeter at 2; = —103 m. 

The Forward Muon Detector (FMD) and the Proton Remnant Tagger (PRT) are sensitive 
to the energy flow in the forward region. They are used to efficiently reject events which do 
not exhibit a rapidity gap between the X system and the proton dissociation system Y. The 
FMD is located at z = 6.5 m and covers a pseudorapidity range of 1.9 < 77 < 3.7. It may also 
detect particles produced at larger 77 due to secondary scattering within the beam pipe. The PRT 
consists of a set of scintillators surrounding the beam pipe at 2 = 26 m and covers the region 
6 < ?7 < 7.5. 

3.2 Event Selection 

The data used in this analysis correspond to an integrated luminosity of 51.5 pb~^ taken in 
the 1999 and 2000 running periods, in which HERA collided protons of 920 GeV energy 
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with positrons of 27.5 GeV. The data are collected using a trigger which requires the scat- 
tered positron to be detected in the SPACAL calorimeter and at least one track of transverse 
momentum above 0.8 GeV to be recorded in the central jet chamber. In the off-line analy- 
sis, the scattered positron is selected as an electromagnetic SPACAL cluster with an energy 
Ee > 8 GeV and polar angle 156° < 9e < 176°. These requirements are well matched to the 
chosen kinematic range of 4 < < 80 GeV^ and 0.1 < y < 0.7. Background from photopro- 
duction, where the positron scatters unobserved at small angles and a particle from the hadronic 
final state is misidentified as the scattered positron, is suppressed by the requirement that the 
difference between the total energy and longitudinal momentum reconstructed in the detector, 
E — pz, must be larger than 35 GeV. Background not related to ep collisions is reduced by 
restricting the z position of the event vertex to lie within 35 cm of the average ep interaction 
point. 

Diffractive events are selected by the absence of hadronic activity above noise threshold in 
the most forward part of the LAr calorimeter (rj > 3.2) and in the FMD and PRT. This selection 
ensures that the rapidity gap between the systems X and Y spans more than four units between 
7] = 3.2 and 7.5. In addition the restriction xp < 0.03 is imposed to limit the contribution from 
secondary reggeon exchanges and to ensure good acceptance. 

The hadronic system X is measured in the LAr and SPACAL calorimeters and the central 
tracking system. Calorimeter cluster energies and track momenta are combined into hadronic 
objects using an algorithm which avoids double counting [32]. Jets are formed from the hadronic 
objects, using the inclusive kx cluster algorithm [25] with a distance parameter of unity in the 
photon-proton rest frame. At least two jets are required with transverse momenta in the 
centre of mass frame of Pxjeti > 5.5 GeV and PT,jet2 > 4 GeV for the leading and sub-leading 
jet, respectively. Asymmetric cuts on the jet transverse momenta are chosen to facilitate com- 
parisons with NLO QCD predictions. The axes of the jets are required to lie within the region 
— 1.0 < 7] jet < 2.0 in the laboratory frame, well within the acceptance of the LAr calorimeter. 
After all cuts 2723 diffractive dijet events are selected. 

3.3 Kinematic Reconstruction 

The energy E^ and polar angle 9^ of the scattered positron are measured using the SPACAL and 
the reconstructed vertex position. The inelasticity y and photon virtuality Q"^ are determined 
according to 



where E^ is the positron beam energy. The energy and momentum of the hadronic system X 
are reconstructed from the observed hadronic objects and the invariant mass Mx is computed 
from this information. The invariant mass of the dijet system is given by 



y 
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with pjeti and pjet2 being the four-momenta of the leading and sub-leading jet, respectively. The 
observables xjp and zp are reconstructed according to 



Xip 
Z]p 



Ml 



ys 

Ml + Q2 



'12 

Mi 



3.4 Monte Carlo Simulations and Fixed Order QCD Predictions 

Monte Carlo simulations are used in the analysis to correct the data for detector effects. For 
events generated with Monte Carlo programs, the HI detector response is simulated in detail 
using GEANT [33] and the events are subjected to the same analysis as the data. Events are 
generated using the RAPGAP program [34] which simulates the process ep eXp, assuming 
proton vertex factorisation (see section |5^ . Leading order matrix elements for the hard QCD 
sub-process are convoluted with DPDFs, taken at the factorisation scale fif = \/p\ + Q^, 
where pt is the transverse momentum of the emerging hard partons relative to the collision axis 
in the 7*j9 centre of mass frame. A preliminary version of the 'HI 2006 DPDF' fit [35] is used 
to simulate pomeron and sub-leading reggeon exchanges. Higher order effects are simulated 
using parton showers [36] in the leading logarithm approximation. The Lund string model 
[37,38] is used for hadronisation. QED radiative corrections are applied using the HERACLES 
program [39]. Processes with a resolved virtual photon are also included, with the structure of 
the photon given by the SAS-2D parameterisation [40]. 

The background due to non-diffractive deep-inelastic scattering is estimated and accounted 
for using the RAPGAP Monte Carlo program in its inclusive mode. The parameters are chosen 
to be similar to the ones used for the generation of the diffractive sample discussed above. The 
inclusive simulation uses the CTEQ5L parton densities of the proton [41]. 

In diffractive DIS measurements using the present technique the system Y does not nec- 
essarily consist only of a single proton, but may also be a low mass dissociative system. 
The DIFFVM program [42] includes a sophisticated treatment of the dissociating proton. It 
is used to study the response of the forward detectors to low mass proton dissociation sys- 
tems (mp < My < 5 GeV). The non-resonant part of the My distribution is modelled with 
da/dMy oc (1/My)^ °^, while the t dependence follows an exponential decrease: da/dt oc e** 
with b = 1.6 GeV~^. This parameterisation is motivated by measurements of diffractive vector 
meson production at HI [43]. Proton dissociation processes with My > 5 GeV are included in 
the treatment of non-diffractive background with RAPGAP as discussed above. 

NLO QCD predictions for the dijet cross sections are calculated at the parton level using the 
NLOJET-I-+ program [44] in slices of xp, assuming proton vertex factorisation. The resulting 
cross sections are converted to the stable hadron level by factors extracted from the RAPGAP 
Monte Carlo mo del in the dif fractive mode. The renormalisation and factorisation scales are set 
to fir = fJ'f = ^Jprjeti + Q^' where j^^-^ is the transverse momentum of the leading jet in the 
7*p centre of mass frame. The NLOJET-i-i- calculation uses parton densities obtained from a 
NLO QCD analysis of inclusive diffractive scattering at HI [5]. That publication provides two 
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4 
0.1 


< Q' < 

< y < 


80 GeV^ 
0.7 






Xjp < 


0.03 


Diffractive Selection 




My < 


1.6 GeV 






\t\ < 


1 GeV^ 


Jet Selection 


-3 


PT,jetl> 
PT,jet2> 
< V]et < 


5.5 GeV 
4 GeV 




Table 1 : The kinematic domain in which the cross sections are measured at the level of jets of 
stable hadrons. The jets are reconstructed using the inclusive kx algorithm as described in the 
text. Variables marked with a are evaluated relative to the collision axis in the 7*p centre of 
mass frame. 



sets of parton densities, HI 2006 DPDF fit A and fit B, which differ in the parameterisation of 
the gluon density. A steeper fall-off in the gluon density at high 2;jp is obtained for fit B than for 
fit A, while the quark densities agree within the uncertainties. Both DPDF sets provide a good 
description of the inclusive diffractive DIS data. 

The experimental and theoretical uncertainties on the DPDFs in [5] are propagated to the 
dijet prediction via an eigenvector decomposition of the error sources according to the method 
presented in [48]. The deviations from the nominal prediction are added in quadrature to obtain 
the uncertainty on the dijet prediction due to DPDF uncertainties. Alternative hadronisation 
corrections are extracted from the POMWIG Monte Carlo model [45], which uses cluster frag- 
mentation [46, 47] to describe hadronisation. The difference between the nominal and alterna- 
tive hadronisation corrections is taken to be the hadronisation uncertainty on the NLO QCD 
prediction. To account for the uncertainty due to the missing higher orders in the calculation, 
the renormalisation and factorisation scales are varied by common factors of 2 and 0.5 with 
respect to the nominal prediction. 

3.5 Cross Section IMeasurement 

The measured differential dijet cross sections are defined at the level of stable hadrons in the 
kinematic region specified in table [TJ A correction of typically 20% is applied to account for 
detector acceptances, inefficiencies and migrations between measurement bins using the RAP- 
GAP 3.1 Monte Carlo program. This simulation gives a reasonable description of the shapes 
of all data distributions. According to the simulation, the detector level observables are found 
to be well correlated with the observables at hadron level. The cross sections are corrected to 
the QED Bom level using the HERACLES interface to the RAPGAP Monte Carlo program. 
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The small background contribution from non-diffractive deep-inelastic scattering is statistically 
subtracted using the Monte Carlo sample introduced above. 

The cross section definition for this study is chosen to include all events with My < 1.6 GeV 
and \t\ < 1 GeV^ as in [2,5, 16-18]. As My and \t\ are not measured directly, the effects of 
migration across these boundaries must be estimated. Migrations from large My > 5 GeV and 
xjp > 0.2 are corrected for using RAPGAP in inclusive mode. Smearing across the My = 
1.6 GeV boundary of events with My < 5 GeV is evaluated with the DIFFYM [42] simulation 
of proton dissociation, following [2]. 

3.6 Systematic Uncertainties 

The systematic uncertainties are evaluated separately for each measurement bin, except for 
uncertainties on global correction factors. The following sources of uncertainty are determined 
to be largely correlated between bins: 

LAr calorimeter energy scale: The energy scale of the LAr calorimeter response to hadrons 
is varied by ±4% in the simulation, which causes a variation of the total cross section by 
'^1% and slightly larger uncertainties in individual measurement bins. 

Track Momenta: The contribution of the track momenta to the X system is varied by ±3%, 
resulting in a total cross section uncertainty of around 3%. 

Luminosity: The measurement of the integrated luminosity has an uncertainty of 1.5%. This 
translates directly into a 1.5% uncertainty on the cross section. 

FMD noise: The cross section is corrected for the fraction of events rejected due to noise in 
the forward muon detector. A global correction factor is determined from a sample of 
randomly triggered events and is found to be (1.2 ± 0.4)%. The uncertainty on this 
correction factor leads to an overall normalisation uncertainty of 0.4%. 

ajj^-migration: The estimated number of non-diffractive background events which migrate 
into the sample from the unmeasured region xp > 0.03 or My > 5 GeV is varied by 
±50%, leading to a total cross section uncertainty of 1%. 

My and |t| migrations: The systematic uncertainties connected to migrations over the My 
and \t\ limits are assessed following the method of [5], giving a total uncertainty of 5%. 

Rapidity gap selection inefficiency: A fraction of the events in the kinematic range specified 
in tableiHgive rise to hadronic activity at pseudorapidities larger than allowed by the r]max 
cut in the LAr calorimeter or in the forward detectors and is thus lost. The correction for 
this effect relies heavily on the RAPGAP simulation to describe the forward energy flow 
of diffractive events. The forward energy flow in diffractive DIS is investigated with a 
sample of elastically scattered protons detected in the forward proton spectrometer of 
the HI detector [49]. The study finds the RAPGAP model to describe these migrations 
to within 30% [50]. The effect of this uncertainty on this measurement is estimated by 
reweighting all events in the signal simulation which do not pass the forward detector cuts 
by ±30%. This translates into an uncertainty of ^5°% on the total cross section. 
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The remaining systematic uncertainties, described below, show significantly less correlation 
and are thus treated as uncorrected between measurement bins. 

Positron energy: The energy of the scattered positron is known to within 2% at i^g = 8 GeV, 
falling linearly to 0.3% at ii^e = 27.5 GeV. This translates into a 2% uncertainty on the 
total cross section. 

Positron angle: The uncertainty in the polar angle 9e of the scattered positron is 1 mrad. This 
contributes an uncertainty of 1% to the total cross section. 

Trigger efficiency: The average difference between the trigger efficiency as extracted from the 
Monte Carlo simulation and from the data using monitor trigger samples is taken as the 
uncertainty on the trigger efficiency, which is around 1%. 

Unfolding uncertainties: To evaluate the model dependence of the correction from the detec- 
tor to the hadron level, key kinematic dependences of the Monte Carlo simulation are 
reweighted within the limits imposed by the present data. The following distributions are 
varied: xp by x^''^, pT by |t| by e^^* '^^^ and y by ?/^°'^. The largest uncer- 

tainty is introduced by the pT reweighting (typically 4%) followed by a;jp (3%), while the 
reweights in the two other variables have rather small effects. 

The largest contributions to the systematic errors on the cross sections arise from the uncer- 
tainty in the LAr calorimeter energy scale, from unfolding uncertainties and from the rapidity 
gap selection inefficiency. The overall uncertainty on the total cross section is The un- 

certainties on individual measurement bins are slightly larger. 

4 Dijet Results 

The integrated cross section in the kinematic range specified in table[T]is determined to be 

a^^^''{ep eXY) = 52 ± 1 (stat.) tl (syst.) pb. 

When this measurement is translated to the kinematic range of the previous HI result [18] 
(i.e. after correcting for the different proton beam energies, y-range and J^^^-ranges), the 
two results are compatible within the uncertainties. The total cross section can be compared to 
the NLO QCD predictions based on the two sets of DPDFs determined from inclusive diffrac- 
tion [5]: 

a^^^^'im 2006 DPDF fit A) = 75 tH (scale unc.) ± 7(DPDF) pb, 
^2jets^^^ 2006 DPDF fit B) = 57 (scale unc.) ± 8(DPDF) pb. 

The scale uncertainty is derived by simultaneously varying i^i j and /if. by common factors of 2 
and 0.5. Whilst both predictions are compatible with the measurement, the central result of fit 
A overestimates the cross section by ~ 40%. 

Differential dijet cross sections are shown in figures [2] to [6] and tabulated in tables l4l to [TOl 
Cross sections as a function of y, xjp, Prjeti ^Vjets = IVjeti ~ Vjet2\ shown in Figured 
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and are compared to NLO QCD predictions. Differential cross sections as a function of zjp axe 
shown in Figure [3l The NLO QCD prediction for the highest bin in zp is not shown due to 
problems in evaluating the hadronisation correction^. 

The prediction based on HI 2006 DPDF fit B describes the shapes of all distributions well, 
whereas some discrepancies are apparent between the fit A and the data. The largest differences 
between the shapes of the two predictions can be seen in 2^ and y, which are correlated through 
the kinematics. The discrepancies between Fit A and the data are most prominent in the region 
of high Z]p {zjp > 0.4), where the prediction is clearly too high. The good agreement in the xjp 
distribution between the dijet data and the predictions indicates that the pomeron flux (which 
governs this distribution) for jet production does not differ significantly from the flux describing 
inclusive diffraction. The shapes of the A?]*^^^ and Prjeti distributions are determined by the 
hard scattering matrix elements and are rather insensitive to the DPDFs. The agreement in 
these distributions shows that the NLO QCD computation, which uses boson gluon fusion as 
the dominant process, is adequate to describe dijet production in this kinematic regime. 

The large difference between the two predictions at high zp reflects the large uncertainty on 
the gluon density in this range as determined from inclusive data alone. Figure [3] also indicates 
the sensitivity of the dijet data to the gluon density at large zjp. To test factorisation in a region 
where the gluon density is well determined from the inclusive data, the dijet cross section is 
also measured in the reduced kinematic domain of zp < 0.4. The results are shown in Fig- 
ured and are compared with predictions based on the HI 2006 DPDF fits. In this kinematic 
region both fits agree well with the dijet data, supporting the notion of QCD factorisation within 
uncertainties. 



5 Combined NLO QCD Fit 

A NLO QCD fit is used to determine the diffractive quark singlet and gluon densities. This 
combined fit uses both the measurements of the diffractive dijet cross sections presented in 
this paper and the measurement of the inclusive diffractive cross section presented in [5]. The 
combined fit shall henceforth be referred to as 'HI 2007 Jets DPDF'. 



5.1 Data Sets 

Assuming the factorisation hypothesis, the differential dijet cross section as a function of zjp is 
used in the fit in four bins of Q'^ +pT,jeti^ which is taken to be the scale variable. These measured 
cross sections are shown in figure [5] and tables |9] and [10] for dijets at the stable hadron level in 
the kinematic range specified in table [T] The fit also includes the measurements of inclusive 
diffraction obtained by HI in [5], which are presented in the form of the reduced diffractive 
deep-inelastic scattering cross section ar^^\ defined through 

'in some cases the Lund string fragmentation algorithm turns the entire system X into just two mesons. This 
leads to events having zjp ~ 1 at the hadron level independently of their parton level zp and to corresponding 
migration problems. 
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where = 1 + (1 — y)^. In leading order the reduced cross section a^'^^^ is identical to 
^D(3) rpj^g small influence of the longitudinal structure function is included here via 

its NLO dependence on the DPDFs. Following the treatment in [5] only data in the range 
g2 > 8.5 GeV^ Mx > 2 GeV and (5 < 0.8 are included in the fit. Figures |7] and [8] show the 
inclusive data points in the form of the product xjp ■ a!!'^^^ [xjp, j3, Q"^). 



5.2 Fit Ansatz 



The DPDFs /ij, xjp, t) are parameterised following the fit procedure of the inclusive 

analysis [5]. They are factorised into a pomeron flux f]p/p{x]p,t) and parton densities of the 
pomeron fi(z, j^j) using the proton vertex factorisation ansatz 



f[^{z,n'j,xip,t) = fip/p{xip,t) ■ fi{z,^i)). 

The parton densities fi are modelled as a singlet distribution /i^) consisting of the three 
light quark and corresponding antiquark distributions, which are all assumed to be of equal 
magnitude, and a gluon distribution g{z, jj^j). Here z is the longitudinal momentum fraction of 
the parton entering the hard subprocess with respect to the diffractive exchange, such that z = 
(3 = x/x]p and z = zjp for the lowest order quark parton model process in inclusive diffraction 
and for dijets, respectively. The parton densities fi(z, j^j) are parameterised at a starting scale 
of /ij Q = 2.5 GeV^ and are evolved to higher factorisation scales using a numerical solution of 
the NLO DGLAP evolution equations. The singlet and gluon distributions are parameterised at 
the starting scale as 



f,iz,fil,)=A-z^^-{l-zf'. 

The parameterisation of the singlet density is thus identical to that used in the analysis of in- 
clusive diffraction [5]. The parameterisation of the gluon density differs in that the HI 2006 
DPDF fit A omits the factor z^3i^°", while fit B omits both z^m^on ^nd (1 — z)^ai^''^\ In the HI 
2007 Jets DPDF fit, where the dijet data additionally constrain the gluon, the of the fit is 
significantly reduced by the inclusion of the factor z^sinon . 

The pomeron flux is parameterised as in [5] using a form motivated by Regge theory: 



\ 2ajp(t)-l 

— e^^*. 

The normalisation parameter Apis defined as in [5]. The pomeron trajectory aip{t) is assumed 
to be linear: 



C(r{t) = ajp(O) + a'jp ■ t. 
13 



Parameter 


Value 


Source 


a'jp 
B]p 
o:]r{0) 

Bm 

rric 

rub 


0.06+°J^ GeV-2 
5.5^°:^ GeV-2 
0.5±0.1 

0.3l°i GeV-2 
1.61?-^ GeV-2 
1.4±0.2 GeV 
4.5± 0.5 GeV 
0.118± 0.002 


[52] 
[52] 
[2] 
[52] 
[52] 
[53] 
[53] 
[53] 



Table 2: Fixed parameters and associated uncertainties used in the HI 2007 Jets DPDF. 



For comparison with the data, all DPDFs are integrated over the measured range |t| < 1 GeV^. 
To properly describe the data, especially at high xjp, it is necessary to include a sub-leading ex- 
change (the so called reggeon, ]R, for details see [5]). This contribution is assumed to factorise 
similarly to the pomeron, so that the definition of the diffractive parton densities is modified to 

ft^iz, /ij, xjp, t) = f]p/p{x]p, t) ■ fi{z, ii)) + rijR ■ f]R/p{xiP, t) ■ f^{z, ii)). 

The reggeon flux f]R/p{x]p, t) is parameterised in the same way as the pomeron flux. The parton 
densities ff^{z, Q"^) are taken from a parameterisation of pion structure function data [51]. The 
free parameters of the fit are the six parameters of the initial parton densities, ajp(O) and the 
normalisation of the reggeon flux rijR. All other parameters are fixed using the same values and 
uncertainties as in [5] as listed in table [2| 

5.3 Fit Procedure 

The fit is performed by minimisation of a function, defined similarly to that in [5]. At each 
step of the minimisation procedure, the predictions for a^^^^ are calculated at NLO in the MS 
renormalisation scheme with the QCDFIT program [54,55]. For the prediction of the dijet cross 
section the combined fit uses the 'matrix method' introduced by the ZEUS collaboration [56], 
together with the NLOJET-i-i- program. This procedure has been shown to yield results which 
agree with direct NLOJET-i-i- predictions in the selected fit range to better than 2% after one 
iteration of the input DPDFs (for details see [57]). Whereas the NLOJET-i-i- calculation employs 
a massless heavy flavour scheme, the prediction for cr^ is performed with massive charm and 
beauty quarks. However, for the dijet data the hard scale is typically much larger than the charm 
mass (nl / > 29 GeV^ ^ m^), so little effect is expected. This is confirmed by performing fits 
to the inclusive data alone in both schemes, resulting in very similar gluon densities in the 
range to which the dijets are sensitive. 

The inclusive and dijet data sets are statistically independent and the correlations between 
the two measurements through the systematic uncertainties are small. The function treats the 
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combined statistical und uncorrelated systematic errors for each data point in the usual way and 
also takes account of correlated uncertainties within the inclusive or dijet data sets by allowing 
variations in the corresponding systematic error sources at the expense of increases in the 
variable [58]. As in [5], there are ten such error sources for inclusive data. The correlated errors 
on the jet cross sections are treated via a single additional parameter in the finction 

Besides the uncertainties related to the cross section measurements, the extracted DPDFs are 
affected by uncertainties in the fit procedure and its theory input. The fit errors include the rather 
small effects of the uncertainties on the input parameters as given in table [2l The uncertainty in 
the relative scale choice between the inclusive and dijet data is estimated by varying the scale for 

the dijet data between 2 ■ ^Jo^^^rp^^^i and 0.5 ■ y'^QM-pp"^ whilst keeping Q as the scale for 
the inclusive data. In addition the effects of changing the fit range m (excluding dijet events 
with Z]p < 0.2) or the starting scale fij^ (using 3.5 GeV^ instead of 2.5 GeV^) are evaluated 
and included in the presented uncertainties. To assess the dependence of the final fitted parton 
densities on the hadronisation correction applied to the dijets, alternative correction factors 
extracted from the POMWIG Monte Carlo model are used. The deviation from the nominal fit 
result is included in the theoretical parton density uncertainties. The largest theoretical error 
contribution to the fitted gluon density at high zjp comes from the uncertainty in the relative 
scale of the two data sets. 



6 HI 2007 Jets DPDF Fit Results 

The fit results for the free parameters are summarized in table |3l The fit describes the data well 
as indicated by the overall value of /ndf = 196/218, which splits into = 27 for the 36 dijet 
data points and = 169 for the 190 a^''^^ data points. Thus the partial for the inclusive 
data is slightly larger in the combined fit than in the fits to cr^*^'^^ from [5], where = 158 
(164) for the HI 2006 DPDF fit A (B), indicating a small remaining tension between the two 
data sets. The parameter Cgiuon, determining the gluon density behaviour at high values of z, 
is positive in the combined fit in accordance with the expectation that the gluon density should 
not be singular for 2; ^ 1. This behaviour is different from the HI DPDF fit A, where Cgiuon 
is determined to be negative and the gluon density is artificially suppressed at very high z using 
an additional exponential factor. 

The dijet cross sections are well described by the predictions based on the HI 2007 Jets 
DPDF as shown in figures [5] and [6l Figures |7] and [8] show the measurements of af^^^^ as a 
function of for different values of /? and xp, together with the NLO predictions based on the 
HI 2007 Jets DPDF fit. The results of fits A and B to the inclusive data alone are also shown. 
A very good description is obtained with all three fits. 

The diffractive gluon distribution and the quark singlet distribution are shown in figure |9]for 
scales of fij = 25 GeV^ and /i^ = 90 GeV^, together with the resuks of fits A and B of the 

stand-alone analysis of a,?^^'' . The error bands indicate the uncertainties due to experimental 
sources and the theoretical errors inherent in the fit procedure. The uncertainties on the quark 
distribution and on the gluon distribution at low Z]p are dominated by the experimental uncer- 
tainties, while the uncertainty on the gluon density at high zjp recieves sizeable contributions 
from both experimental and theoretical sources. 
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Parameter 


Fit Value 
(HI 2007 Jets DPDF) 


aip{0) 


1.104 ±0.007 




1.3xl0-3±0.4 X 10-3 


-^gluon 


0.88 ±0.17 


Bgluon 


0.33 ±0.10 


Cgluon 


0.91 ±0.18 


-^quark 


0.13 ±0.02 


B quark 


1.5 ±0.12 


C quark 


0.51 ±0.08 




196/218 



Table 3: HI 2007 Jets DPDF fit parameters obtained from the combined fit to the diffractive 
inclusive and dijet data. Only the experimental uncertainties are given. 



The combined fit constrains both the diffractive gluon and quark densities well and for the 
first time with comparable precision in the complete range 0.05 < zjp < 0.9. At high zjp the 
resulting gluon density differs significantly from that of HI 2006 DPDF fit A, but is compatible 
with fit B [5]. Good agreement is seen between all three fits for the singlet quark density and 
the gluon density at low ^jp. The values of ajp(O) = 1.104±0.007 and rijR = (1.3±0.4) x IQ-^ 
are compatible within experimental uncertainties with the value extracted in HI 2006 fit B. The 
uncertainties on these parameters are not significantly decreased by the inclusion of the dijet 
data compared to the determination from the inclusive data alone. 

In figure[lO]the DPDFs as determined by HI are compared with the results of an independent 
analysis [15], where parton densities are derived from the same inclusive diffractive data [5]. 
A hybrid theoretical framework is used which combines aspects of collinear factorisation and 
a perturbative two-gluon-exchange model [59-62]. Most of the dijets events are produced via 
BGF-type processes as in figure \T\ At high P, there is an additional contribution in which the 
perturbative two-gluon state participates directly in the hard interaction via photon-pomeron 
fusion, leading to a modified evolution equation for the DPDFs. The resulting DPDFs agree 
reasonably well with the HI 2007 Jets DPDF and with the HI 2006 DPDF fit B. 

Measurements of diffractive charm production by HI have also been compared to predic- 
tions based on the DPDFs presented in this paper [63]. Whilst overall good agreement is 
obtained, the statistical accuracy of the charm measurement limits its power to discriminate 
between different DPDF sets. 



7 Conclusion 

Cross sections for dijet production in diffractive deep-inelastic scattering are measured with 
improved precision compared to earlier analyses. Single and double differential cross sections 
are presented in a variety of variables sensitive to the underlying dynamics of hard diffraction. 
NLO QCD predictions based on diffractive parton densities extracted from measurements of in- 
clusive diffractive deep-inelastic scattering describe the data well in the kinematic region where 
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the gluon density is reliably constrained by the inclusive measurements. This agreement con- 
firms the validity of QCD factorisation and thus the applicability of diffractive parton densities 
evolving according to the DGLAP equations. 

A combined fit to diffractive inclusive and dijet data is performed, using NLO QCD calcu- 
lations based on QCD factorisation and DGLAP evolution. Both data sets are described well 
by the fit. The inclusion of the dijet data allows the simultaneous determination of both the 
diffractive gluon and the singlet quark distribution with good and comparable accuracy in the 
range 0.05 < zjp < 0.9. This is the first reliable determination of the diffractive gluon density 
up to large momentum fractions. 
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Figure 2: Cross sections for diffractive dijets, differential in y, \ogxjp, Prjeti ^^'^ ^Vjets com- 
pared to NLO predictions based on the parton-densities from the HI 2006 DPDF fits [5]. The 
data are shown as black points with the inner and outer error bars denoting the statistical and 
quadratically added uncorrelated systematic uncertainties, respectively. The hatched band in- 
dicates the correlated systematic uncertainty. The dashed line shows the NLO QCD prediction 
based on the HI 2006 DPDF fit B, which is surrounded by a dark shaded band indicating the 
parton density and hadronisation uncertainties. In the light shaded band the scale uncertainty 
is added quadratically to the parton density and hadronisation uncertainties. The dotted line 
represents the NLO QCD prediction based on the HI 2006 DPDF fit A. 
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Figure 3: Cross section for diffractive dijets, differential in zp compared to NLO predictions 
based on the parton-densities from the HI 2006 DPDF fits [5]. The data are shown as black 
points with the inner and outer error bars denoting the statistical and quadratically added uncor- 
rected systematic uncertainties, respectively. The hatched band indicates the correlated system- 
atic uncertainty. In the left panel the data are compared to the NLO QCD prediction based on 
the HI 2006 DPDF fit A (dotted line) and in the right panerl to the prediction based on the HI 
2006 DPDF fit B (dashed line). The lines are surrounded by a dark shaded band indicating the 
parton density and hadronisation uncertainties. In the light shaded band the scale uncertainty 
is added quadratically to the parton density and hadronisation uncertainties. The prediction for 
zp > 0.9 is not shown since the hadronisation corrections for this bin cannot be determined 
reliably. 
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Figure 4: Cross sections for diffractive dijets restricted to zp < 0.4, differential in y, \ogxp, 
Prjeti ^Vjets compared to NLO predictions based on the parton-densities from the HI 2006 
DPDF fits [5]. The data are shown as black points with the inner and outer error bars denoting 
the statistical and quadratically added uncorrelated systematic uncertainties, respectively. The 
hatched band indicates the correlated systematic uncertainty. The dashed line shows the NLO 
QCD prediction based on the HI 2006 DPDF fit B, which is surrounded by a dark shaded band 
indicating the parton density and hadronisation uncertainties. In the light shaded band the scale 
uncertainty is added quadratically to the parton density and hadronisation uncertainties. The 
dotted line represents the NLO QCD prediction based on the HI 2006 DPDF fit A. 
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Figure 5: Cross section for diffractive dijet production doubly differential in zp and the scale 
+ PT,jeti- The data are shown as black points with the inner and outer error bars denoting 
the statistical and quadratically added uncorrected systematic uncertainties, respectively. The 
hatched band indicates the correlated systematic uncertainty. The solid line shows the NLO 
QCD prediction based on the HI 2007 Jets DPDF. Data points in the highest Z]p bin were not 
included in the fit since the hadronisation corrections cannot be evaluated reliably. 
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Figure 6: Cross sections for diffractive dijet production differential in the variables y, \ogxp, 
Prjeti ^Vjets- The data are shown as black points with the inner and outer error bars denot- 
ing the statistical and quadratically added uncorrected systematic uncertainties, respectively. 
The hatched band indicates the correlated systematic uncertainty. The solid line surrounded 
by the shaded band shows the NLO QCD prediction based on the HI 2007 Jets DPDF, where 
the band denotes the scale uncertainty derived by varying the renormalisation and factorisation 

scale ^ = JQ"^ + PTjeti by factors of 2 and 0.5. 
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Figure 7: The dependence of the diffractive reduced cross section a!^^^^ muhiplied by at Xip=0.001 (left) and xjp=0.003 (right) at 
various values of p. The cross sections are multiplied by powers of 3 for better visibility. The data points are taken from the publication [5]. 
The inner and outer error bars on the data points represent the statistical and total uncertainties, respectively. Only data points included in 
the DPDF fits are shown. The data are compared to NLO QCD predictions based on the HI 2007 Jets DPDF, which are shown as solid lines. 
The dashed and dotted lines indicate the predictions of the HI 2006 DPDF fit A and B, respectively. 
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Figure 8: The Q'^ dependence of the diffractive reduced cross section cr,?''^^ multiplied by Xjp at X]p=O.Ol (left) and xjp=0.03 (right) at 
various values of p. See caption of figure|7]for further details. 
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Figure 9: The diffractive quark density (top) and the diffractive gluon density (bottom) for two 
values of the squared factorisation scale /i^: 25 GeV^ (left) and 90 GeV^ (right). The solid line 
indicates the HI 2007 Jets DPDF, surrounded by the experimental uncertainty (dark shaded 
band) and the experimental and theoretical uncertainties added in quadrature (light shaded 
band). The dotted and dashed lines show the parton densities corresponding to the HI 2006 
fit A and fit B from [5], respectively. 
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Figure 10: The diffractive quark density (left) and the dilfractive gluon density (right) for the 
factorisation scale f-ij = 25 GeV^ . The solid line indicates the HI 2007 Jets DPDF. The dotted 
and dashed lines show the parton densities corresponding to the HI 2006 fit A and fit B from [5], 
respectively. The dashed-dotted line shows the DPDFs as determined by Martin, Ryskin and 
Watt in [15]. 
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Table 4: Bin averaged differential cross sections of diffractive dijet production at the hadron 
level (corrected to the QED Born level) and the corresponding uncertainties as a function of xjp. 
The corrections applied to the NLO prediction for hadronisation and the associated uncertainty 
are also given. 
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Table 5: Bin averaged differential cross sections of diffractive dijet production at the hadron 
level (corrected to the QED Bom level) and the corresponding uncertainties as a function of y. 
The corrections applied to the NLO prediction for hadronisation and the associated uncertainty 
are also given. 
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Table 6: Bin averaged differential cross sections of diffractive dijet production at the hadron 
level (corrected to the QED Born level) and the corresponding uncertainties as a function of ^jp. 
The corrections applied to the NLO prediction for hadronisation and the associated uncertainty 
are also given. No hadronisation correction is given for the highest zjp bin since it cannot be 
evaluated reliably. 
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Table 7: Bin averaged differential cross sections of diffractive dijet production at the hadron 
level (corrected to the QED Born level) and the corresponding uncertainties as a function of 
Prjeti- The corrections applied to the NLO prediction for hadronisation and the associated 
uncertainty are also given. 
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Table 8: Bin averaged differential cross sections of diffractive dijet production at the hadron 
level (corrected to the QED Born level) and the corresponding uncertainties as a function of 
^Vjets- The corrections applied to the NLO prediction for hadronisation and the associated 
uncertainty are also given. 



31 









29GeV2 < Q^+P^,,,,i 


< 50 GeV^ 










zp 




d^a/dzpdfi^ 


6tot. 




^uncorr. 






hadr. 








[pb/GcV^] 


[pb/GeV^] 


[pb/GcV^] 


[pb/GcV^] 


[pb/GcV^] 




corr. 




0.0 - 


0.1 


0.16 


0.10 


0.05 


0.04 


0.07 


1.32 


± 


0.04 


0.1 - 


0.2 


1.10 


0.30 


0.15 


0.18 


0.19 


0.99 


± 


0.35 


0.2 - 


0.3 


1.24 


0.25 


0.17 


0.10 


0.15 


1.09 


± 


0.11 


0.3 - 


0.4 


1.16 


0.24 


0.17 


0.07 


0.15 


0.97 


± 


0.14 


0.4 - 


0.5 


1.12 


0.23 


0.18 


0.06 


0.13 


1.08 


± 


0.01 


0.5 - 


0.6 


0.61 


0.17 


0.10 


0.07 


0.11 


1.11 


± 


0.10 


0.6 - 


0.7 


0.45 


0.12 


0.09 


0.05 


0.06 


0.91 


± 


0.01 


0.7 - 


0.8 


0.197 


0.071 


0.056 


0.031 


0.030 


0.86 


± 


0.60 


0.8 - 


0.9 


0.042 


0.036 


0.022 


0.024 


0.015 


0.98 


± 


0.50 


0.9 - 


1.0 


0.22 


0.22 


0.18 


0.09 


0.09 




- 










50GeV2 < 


< 70 GeV^ 










zp 




d^a/dzpdii^ 


Stot. 


<5stat. 


^^uncorr. 






hadr. 








[pb/GcV^] 


[pb/GcV^] 


[pb/GcV^] 


[pb/GcV^] 


[pb/GcV^] 




corr. 




0.0 - 


0.1 


0.124 


0.059 


0.047 


0.018 


0.030 


1.21 


± 


0.70 


0.1 - 


0.2 


1.52 


0.32 


0.18 


0.11 


0.25 


1.10 


± 


0.10 


0.2 - 


0.3 


1.91 


0.34 


0.20 


0.15 


0.23 


1.08 


± 


0.12 


0.3 - 


0.4 


1.54 


0.33 


0.17 


0.14 


0.25 


1.14 


± 


0.02 


0.4 - 


0.5 


1.18 


0.23 


0.15 


0.09 


0.15 


1.07 


± 


0.20 


0.5 - 


0.6 


1.09 


0.23 


0.15 


0.08 


0.15 


1.08 


± 


0.12 


0.6 - 


0.7 


0.74 


0.15 


0.12 


0.05 


0.07 


1.13 


± 


0.08 


0.7 - 


0.8 


0.35 


0.12 


0.08 


0.05 


0.06 


0.92 


± 


0.36 


0.8 - 


0.9 


0.20 


0.13 


0.09 


0.07 


0.06 


0.76 


± 


0.65 


0.9 - 


1.0 


0.11 


0.13 


0.07 


0.02 


0.10 









Table 9: Bin averaged double differential cross sections of diffractive dijet production at the 
hadron level (corrected to the QED Bom level) and the corresponding uncertainties as a function 
of zjp in different bins of jj? = Q'^ + Prjeti- The corrections applied to the NLO prediction 
for hadronisation and the associated uncertainty are also given. No hadronisation correction is 
given for the highest zp bin since it cannot be evaluated reliably. 
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Table 10: Bin averaged double differential cross sections of diffractive dijet production at the 
hadron level (corrected to the QED Bom level) and the corresponding uncertainties as a function 
of zjp in different bins of jj? = Q"^ + Pr jeti- The corrections applied to the NLO prediction 
for hadronisation and the associated uncertainty are also given. No hadronisation correction is 
given for the highest zjp bin since it cannot be evaluated reliably. 
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